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A  MOIEL  FOR  WIND  FLOW  IM  AM  IDBALIZED  VEQBTATIVE  CAMOPY 


By 

R.  M.  Cloneoy  W.  D.  Ohastada,  and  J.  7.  Applaby 


OBJECTIVE 


Tha  objaotiva  of  Toric  1-A»0«.li001-B»021»08.  "Mieromtaorology 
(ERDAA),"  !■  to  oonduot  stadlos  daallnc  with  tha  phyaleal  prooaaaaa 
involvad  In  tha  axohanga  of  anargy  batwaan  tha  atnoaphara  and  tha 
aarth*8  aarfaea.  Through  aueh  baaie  raaaaroh,  Inoraaaad  knowladga  of 
ataoapharlo  prooaaaaa  will  raaalt  which  can  ba  aaad  by  appliad  raaaaroh 
aotivitiaa  to  iaprora  tha  daaign  and  taat  eritaria  for  Amy  waapona  and 
ayataaa  and  to  ii^rova  waathar  obaarving  and  foreeaating  for  Amy 
taetioal  opamtiona* 
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Th*  piupoM  of  tho  rtpori  i«  to  prosont  tho  rosulta  of  •  •to^jr 
whioh  ox^nos  tho  ohoraetoristlo  poroMtoro  of  wind  flow  In  on 
Idoollsod  TogototlTo  oonopjr.  Equations  oharaotorising  tho  wind  flow 
in  and  abowo  tho  wogotativo  oanopjr  aro  dorivod  and  oooparod  with  data 
eoUootod  in  and  abovo  varions  eropa. 

Obsonrationa  takon  in  a  corn  fiold  indioato  a  aarkod  aladlaritgr 
to  tho  flow  doriwod  for  the  idoal  eanopj.  ilthoogh  tho  atili^  of 
tho  approaoh  was  not  prowon*  tho  ooneapt  dooa  aoaw  to  show  prowiao. 

It  is  hopod  this  roport  will  stbaolato  thooght  and  inTOstigation 
which  will  oltlauitolx  load  to  tho  objootiTO  of  oaproaaing  tho  aoro- 
dsmawio  ronghnoas  offsets  in  torws  of  wagotatiwo  eharaotoriatioa. 
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RESEARCH  STDDT 

A  MODEL  FOR  WIND  FLOW  IN  AN  IDEALIZED  VEGETATIVE  CANOPT 
DA  T>dc  l^A-O-llOOl^B-021-08 


INTRODUCTIOW 

To  a  large  extent,  BloroMteorolog7  Is  a  solenee  deroted 
to  the  evaluation  of  physical  processes  at  the  boundary  or  inter¬ 
face  between  the  earth  and  the  atnosphere.  With  regard  to 
atmospheric  turbulence  within  the  boundary  layer,  the  most 
important  single  characteristic  of  the  boundaiy  is  its  effective 
aerodynamic  roughness.  Over  land  surfaces,  with  few  exceptions, 
the  boundary  is  considered  to  be  fully  rough;  that  is,  a  laminar 
sublayer  is  essentially  nonexistent.  Under  natural  conditions, 
the  roughness  over  land  is  primarily  associated  with  vegetation, 
although,  in  urban  areas  man-made  structures  may  dominate. 

For  the  general  application  of  micrometeorologioal  prediction 
models,  it  would  be  most  advantageous  if  the  aerodynamic  roughness 
effect  of  the  boundary  could  be  expressed  in  terms  of  the  height, 
density  and  drag  characteristics  of  the  roughness  elements.  To 
date,  no  such  goal  has  been  achieved.  Consequently,  it  was  the 
objective  of  the  study  reported  here  to  investigate  the  turbulent 
transfer  of  momentum  within  a  vegetative  canopy  and  to  atteqst  to 
develop  a  preliminary  model  which  expressed  the  aerodynamic  rough¬ 
ness  effects  in  the  terms  of  characteristics  of  the  vegetation. 

BACKGROUND 

For  practical  reasons,  it  has  been  the  goal  of  the  majority 
of  workers  in  turbulence  research  to  express  the  turbulent  trans¬ 
fers  in  terms  of  mean  quantities  of  the  velocity  and  its  derivatives 
rather  than  the  variances  and  covariances  of  the  tnibulent  fluctu¬ 
ations  which  are  known  as  the  Reynolds  stresses.  Although  the 
latter  are  fundamentally  sound,  they  consist,  in  essence,  of  only 
statistical  quantities  which  describe  the  fluctuation  kinematics. 

It  is  the  intent  here  to  continue  the  practice  of  e:q>ressing  the 
dymnic  characteristics  of  the  turbulent  transfer  in  tenu  of  the 
mean  quantities. 

Many  of  the  developments  of  models  for  turbulent  transfer 
have  come  about  through  reasoning  on  dimensional  grounds  and  lean 
heavily  on  empirical  quantities  and  relationships.  There  has  been 
only  limited  success  in  pursuing  the  problem  in  this  manner  and  it 
is  safe  to  say  that  a  general  model  for  atmospheric  turbulence 
remains  remote.  On  the  positive  side,  however,  remarkable  success 
has  been  achieved  in  special  cases  with  very  simple  models.  Fore- 
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most  In  this  srss  Is  the  mlxlng«length  hypothssis  sttrlbutsd  to 
Prsndtl.  This  oonoept  viljl  be  used  later  in  the  dlsoussion  of  the 
turtnilent  transfer  within  Ithe  vegetative  canopy,  so  we  shall 
review  briefly  the  eharacterlstloe  of  the  mixing-length  hypothesis 
for  fully  rough,  steady  flow  at  very  high  Reynolds  numbers.  The 
basic  relationship  of  this  hypothesis  Is  as  follows: 


where  t  Is  the  shearing  stress,  p  Is  the  density,  u^^  Is  the  fric¬ 
tion  velocity,  I  Is  the  mixing  length,  u  Is  the  mean  velocity 
parallel  to  the  boundary,  and  a  Is  the  space  coordinate  normal  to 
the  boundary.  If  the  analogy  with  laminar  flow  Is  Invoked  such 
that  the  shearing  stress  Is  proportional  to  the  shear,  we  obtain: 


(2) 


where  K  Is  the  eddy  viscosity. 

The  behavior  of  the  mixing  length  In  the  vicinity  of  a  rough 
boundary  can  best  be  demonstrated  by  citing  observations  of  the 
boundary  layer  of  a  flat  plate  at  very  high  Reynolds  numbers.  These 
observations  show  that,  near  the  boundary,  the  friction  velocity  Is 
essentially  constant  and  the  shear  Is  Inversely  proportional  to 
the  height  above  the  boundary.  Consequently,  the  mixing  length  must 
be  proportional  to  the  height  above  the  boundary.  This  Is  expressed 
as  follows: 


i  -  ka  (3) 

where  k  Is  von  Karman*s  constant.  The  real  success  of  the  mixing- 
length  hypothesis  lies  In  the  fact  that  k  Is,  Indeed,  a  constant 
which  Is  Independent  of  the  scale  of  flow;  that  Is,  equation  (3)  Is 
equally  applicable  to  the  atmospheric  surface  boundary  layer  as 
well  as  liquid  flow  In  pipes. 

Combining  equation  (1)  with  equation  (3)  and  Integrating  re¬ 
sults  In  the  following  velocity  profile  relationship: 

u  ,  1  fin  1  +  c  1 

^  k  [  J  (4) 

where  c  Is  a  constant  of  Integration.  It  Is  through  this  constant 
of  Integration  that  the  roughness  enters  the  equations  of  turbulent 
transfer.  The  most  oonprehenslve  quantitative  expression  for  c  Is 
due  to  Nlkuradse,  (Reference  1).  He  conducted  extensive  experiments 
on  flow  through  pipes  roughened  by  tightly  packed  sand  grains  of 
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eh«raet«ri8tlo  sise  kg.  For  folly  rough  flow,  hl«  results  showed 
that  e  -  In  (30/kg), 

The  expression  "fully  rough  flow"  requires  explanation.  For 
pipe  flow,  it  is  useful  to  define  a  drag  ooeffioient  (Cq)  as 
follows: 

<  -T./p  *lCnB* 

t  (5) 

where  Tq  is  the  shearing  stress  at  the  boundary  and  U  is  the  mean 
flow  velocity  through  the  pipe.  A  fully  ngugh  flow  is  one  in  which 
Cj)  is  independent  of  the  Remolds  nunber,  Ud/v,  where  d  is  the 
diameter  of  the  pipe  and  v  is  the  kinematic  viscosity.  The  fully 
rough  regime  occurs  for  Reynolds  numbers  exoeeding  a  certain  criti¬ 
cal  value.  This  critical  Reynolds  number  increases  for  decreasing 
relative  roughness  (kg/d). 

The  e:q;>eriment8  of  Nikuradse  incorporated  tightly  packed  sand 
grains  and  no  attenpt  was  made  to  vary  the  density  of  the  roughness 
elements.  Sohlichting  (Reference  2)  conducted  experiments  analo¬ 
gous  to  those  of  Nikuradse  but  with  varying  densities  of  regularly 
shaped  roughness  elements  (spheres  and  hemispheres).  He  expressed 
his  results  in  terms  of  the  effective  sand  grain  roughness  of 
Nikuradse.  As  would  be  anticipated,  the  effective  sand  grain 
roughness  increased  with  increasing  density  of  roughness  elements. 
One  interesting  exception  to  this  rule  was  that,  for  the  most 
dense  arrangement  of  spheres,  the  effective  sand  grain  roughness 
was  less  than  for  the  same  spheres  with  less  packing.  This  result 
can  be  interpreted  as  meaning  the  ^heres  were  so  tightly  packed 
that  the  fluid  "saw"  a  new  *ero  level  displaced  inward  from  the 
actual  wall  of  the  pipe;  that  is,  the  fluid  recognised  only  a  por¬ 
tion  of  the  actual  height  of  the  ^>heres  in  terms  of  roughness. 

In  the  past,  the  treatment  of  roughness  effects  on  the 
atmospheric  boundary  layer  have  largely  followed  the  precedents 
set  in  the  above  consideration  of  fully  rough  pipe  flow  but  with 
some  change  of  parameters.  The  most  general  expression*  for  the 
wind  profile  equation  is  as  follows: 


Equation  (6)  is  applicable  only  to  regions  inasdiately  above  the 
boundaiy  where  the  shearing  stress  is  constant  with  height  and 
effects  of  buoyancy  in  diabatio  conditions  are  negligible.  This 
condition  is  generally  satisfied  in  the  first  two  or  three  meters 
above  the  top  of  the  roughness  elements. 


6 


idMM  D  is  ths  Mro-plsM  dlsplaoMMnt  and  tg  is  ths  rooghnass 
length. 

The  behavior  of  the  sero-plane  displacement  and  the  roughness 
length  is  well  demonstrated  bjr  the  results  of  a  unique  set  of 
experiments  reported  by  Kutsbaoh  (Reference  3)  •  These  experiments 
mere  scswmhat  analogous  to  those  of  Schlichting  described  above. 

The  friction  velocity*  aero-plane  di-splacement,  and  roughnees  length 
mere  evaluated  for  an  enseaible  of  bushel  baskets  of  varying  densl.ty 
patterns  placed  on  an  ice-covered  lake.  The  results  shomed  that 
the  roughness  length  increased  approximately  linearly  mlth  inoress- 
ing  density  of  baskets*  mhile  the  aero-plane  displacement  increased 
mith  the  basket  density  to  the  0.3  pomer.  Although  there  mas 
considerable  scatter  in  the  aero'^lane  displacement  i*elationshlp » 
the  results  indicate  a  rather  orderly  variation  of  the  roughness 
parameters  mith  the  density  of  baskets.  In  another  eoqperiment* 
the  density  of  the  bushel  baskets  mas  maintained  constant*  but  tkM 
height  of  the  obstacles  mas  varied  by  using  successively  one  basket 
and  then  doubling  the  height  by  placing  one  basket  on  top  of  ano-bher. 
As  a  result*  the  zero-plane  displacement  Increased  by  a  factor  oxf 
2.5.  The  difference  from  a  factor  of  2  Is  likely  due  to  a  change 
of  shape  of  the  roughness  element  from  one  Inverted  basket  to  that 
of  an  Inverted  basket  placed  on  top  of  an  up-rlght  basket.  In  a 
qualitative  sense  this  result  Is  compatible  mith  the  results  of 
Schlichting  mhlch  Indicated  that  spheres  are  rougher  than  hemi¬ 
spheres. 

The  bushel  basket  results  Indicate  an  orderly  variation  of 
the  roughness  parameters  mith  the  density  and  height  of  rigid 
roughness  elements.  However,  the  results  are  considerably  less 
definitive  when  the  ensemble  of  roughness  elements  Is  a  uniform 
field  of  vegetation.  This  Is  clearly  Indicated  In  results  publish¬ 
ed  by  Tan  and  Ling  (Reference  4)  and  Stoller  and  Lemon  (Reference  5) 
for  wind  profiles  over  fields  of  wheat,  alfalfa,  and  com.  Briefly, 
their  results  with  com  show  that,  for  a  given  crop  density  and 
height,  the  zero-plane  displacement  decreases  and  the  roughness 
length  Increases  with  Increasing  wind  speed;  in  alfalfa  the  oppo¬ 
site  is  true  with  the  unexpected  result  that  the  shearing  stress 
is  essentially  Independent  of  wind  speed.  It  is  not  difficult  to 
Imagine  the  cause  of  these  changes.  Alfalfa  Is  limber  so  with 
Increasing  wind  speed  the  stems  bend  over  and  the  leaves  tend  to 
orient  themselves  along  lines  of  least  resistance.  The  laying 
over  of  the  plants  partially  seal  off  the  canopy  from  penetration 
of  the  mainstream  above.  The  result  Is  decreasing  roughness  and 
Increasing  zero-plane  displacement  with  Increasing  wind  ^eed  fox* 
a  liidMr  vegetative  cancpy.  Com  Is  characterized  by  seml-rlgld 
stalks  and  limber  but  resilient  leaves.  Increasing  wind  speed 
results  in  greater  penetration  into  the  somewhat  open  canopy  and 
a  flapping  of  the  leaves.  The  result  is  decreasing  zero-plane 
displacement  and  Increasing  roughness  with  Increase  In  wind  speed  in 
com.  So  long  as  the  roughness  length  and  the  zero-plane  dlsplace- 
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nent  are  the  sole  paraaetere  ased*  It  would  appear  that  little 
nore  than  the  qualitative  reasoning  applied  above  eould  be  invoked 
to  characterise  the  roughness  effects  of  a  vegetative  canopy*  The 
point  is  that  the  roughness  length  and  the  aeroplane  displacement 
are  merely  convenient  parameters  to  specify  the  boundary  condition 
of  the  logarithmic  vivA  profile  and  are  only  vaguely  related  to  the 
actual  physical  processes  involved  in  the  roughness  effects*  For 
this  reason  an  entirely  new  approach  was  utilised  in  this  study 
to  evaluate  the  roughness  effects  of  a  vegetative  canopy*  The 
intent  was  to  develop  a  model  for  turbulent  transfer  within  the 
canopy  which  expressed  the  boundary  conditions  for  the  logarithmic 
wind  profile  directly  in  terms  of  the  height,  density  and  drag 
characteristics  of  the  vegetation. 


thb  ideal  CANOFT 


For  preliminary  considerations  of  turbulence  within  the  vegeta¬ 
tion,  let  us  consider  the  steady-state  adiabatic  tux^milent  energy 
equation  (see  Townsend  Reference  6)  in  the  following  fom: 


• ,  j  ®  1«2,3 


u.» 

I 


u» 

j 


axj 


e  =  0 


(7) 


where  the  subscripts  designate  the  three  components  in  a  rectangular 
coordinate  system,  u!.^  represents  the  turbulent  velocity  component 
in  the  i  j  direction,'^  a  i  i  is  the  local  mean  velocity  component  in 
the  i  j  direction,  e  is  tl»  local  dissipation  and  a  bar  designates  a 
time  average*  Equation  (7)  is  derived  from  the  Navier-Stokes  equation 
and  is  assumed  to  apply  to  any  small  localized  region  of  the  fluid 
within  the  vegetative  canopy. 


To  apply  equation  (7)  in  a  practical  way,  let  us  define  a 
volume  V  which  is  a  horizontal  slice  containing  a  representative 
quantity  and  arrangement  of  leaves  and  stalks.  For  this  volume,  we 
can  define  a  gross  mean  velocity  U  by  integration  over  the  volume  of 
the  localised  mean  velocity  components  U|  *  It  is  reasonable  to 
assume  that  U  is  horizontal  and  parallel  with  the  airibient  mean  veloc¬ 
ity  above  the  vegetative  canopy*  Furthermore,  we  define  a  length 
scale  L  which  for  now  will  be  assumed  to  characterise  the  average 
leaf -stalk  arrangement  within  the  volume  V*  Next,  we  non-dimension- 
^ize  the  terms  of  equation  (7)  as  follows*  v',  j  ■  u*j  :  M 
^ii  *  ^ii  7  -  *  ./b^th  the  result  that  equation  (7) 

becomes**  " 


£ 

L 


2lj 


(8) 
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Integrating  equation  (8)  over  the  volume  V,  we  obtain  the  folloirlng 
expression  for  the  mean  dissipation! 

e  »  (c^\  U» 

\  L  /  (9) 


where  the  dissipation  coefficient  C**  of  the  vegetation  within  the 
volume  V  Is  defined  as: 


C* 


sf]! 


dxdjrds 


(10) 


To  complete  the  picture,  we  obtain  the  average  values  of  the  Reynolds 
stresses  within  the  volume  V  by  like  Integration.  If  we  essome,  as 
usual,  that  the  horizontal  derivatives  of  the  gross  mean  quantities 
of  the  volume  V  are  negligible  and  remembering  that  u**  ■  -  u^w*, 
then  equation  (7)  can  be  written  to  describe  the  mean  dissipation  of 
the  volume  V  as  follows: 


u**  M  -  i"  »  (C*  /L)  U» 

h*  (11) 

Rearranging  terms,  we  obtain  the  following  equation: 

Min  U)  C* 

'  L  (12) 

Equation  (12)  represents  a  preliminary  expression  of  turbulent 
transfer  of  momentum  within  a  ve^tatlve  canopy*  In  Its  derivation, 
we  have  Introduced  two  terms— ithe  dissipation  coefficient  C*,  and  the 
scale  length  L.  The  usefulness  of  equation  (12)  depends  upon  whether 
these  two  factors  can  be  evaluated  experimentally  and  expressed  egBlri- 
cally  as  unique  functions  of  the  density  and  structure  (or  species)  of 
vegetation. 


A  dissipation  coefficient  can  be  defined  for  flow  through  a 
fully  rough  pipe  In  the  sane  manner  as  was  done  for  the  drag  coeffi¬ 
cient.  In  the  ease  of  pipes,  the  scale  length  L  Is  the  diameter  of 
the  pipe*  For  sufficiently  high  Reynolds  numbers,  the  dissipation 
coefficient  is  Independent  of  Reynolds  number,  but  Its  magnitude 
depends  vpcn  the  relative  roughness  of  the  pipe  (kg/d).  It  should  be 
noted  that  equation  (9),  which  defines  the  dissipation  coefficient, 
differs  from  equation  (5)  which  defined  the  drag  coefficient  by  inclu¬ 
sion  of  the  scale  length  L.  Thus,  at  very  high  Reynolds  numbers, 
greater  dissipation  occurs  In  small  pipes  than  In  larger  pipes,  even 
though  the  mean  velocities  and  shearing  stresses  nay  be  the  sane*  If 
we  apply  this  reasoning  to  a  vegetative  canopy,  we  would  assume  that 
the  scide  length  in  some  way  measures  the  average  distance  between 
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adjacent  parts  of  the  plants.  If  this  is  the  easey  we  would  expect 
a  higher  degree  of  local  dissipation  associated  with  closely  packed 
plants  such  as  grass  or  alfalfa  as  opposed  to  more  loosely  arranged 
canopies  such  as  a  com  field  or  forest. 

The  terra  (u*/U)  is  somewhat  analogous  to  the  intensity  of  the 
turbulence.  It  is  not  unreasonable  to  assume  that  the  turbulence 
intensity  is  closely  correlated  with  the  drag  coefficient  (see 
equation  (5))  and  thus  the  dissipation  coefficient. 

For  an  ideal  canopy,  let  us  assume  that  the  dissipation  coeffi¬ 
cient  and  the  scale  length  L  are  constant  with  height  in  the  canopy. 
This  is  equivalent  to  assuming  that  the  density,  size,  and  arrange¬ 
ment  of  the  leaves  and  stalks  are  "essentially"  invariant  within 
the  canopy.  Because  of  the  inferred  proportionality  of  the  turbulence 
intensity  to  dissipation  coefficient,  we  further  characterize  the 
ideal  canopy  as  having  a  constant  turbulence  intensity  as  expressed 
by  (u*/U).  Consequently,  the  term  a(ln  U  )/^z  in  equation  (12) 
must  also  be  a  constant.  Thus,  the  velocity  profile  of  the  ideal 
canopy  must  be  an  exponential  function  of  height. 

The  effect  of  roughness  on  the  surface  boundary  layer  becomes 
evident  only  in  the  integrated  equation  which  defines  the  velocity 
profile.  The  profile  equation  is  derived  by  means  of  the  mixing- 
length  hypothesis.  Equation  (12)  was  derived  independent  of  the 
mixing-length  hypothesis,  and  therefore  is  somewhat  incompatible  with 
the  more  common  approach  to  turbulent  transfer. 


Ordway,  Ritter,  Spence,  and  Tan  (Reference  7)  have  suggested 
that  the  turbulent  transfer  within  the  vegetative  canopy  can  be 
pursued  along  the  same  lines  as  the  surface  boundary  layer  provided 
allowance  is  made  for  the  loss  of  momentum  to  the  leaves  and  stalks. 
On  heuristic  grounds,  they  proposed  that  this  momentum  loss  was 
proportional  to  the  square  of  the  local  mean  velocity.  This  is 
analogous  to  equation  (5)  which  defined  the  drag  coefficient  for 
fully  rough  flow  through  a  pipe.  Further,  for  steady  state  condi¬ 
tions  with  no  advective  terms,  the  local  loss  of  momentum  must 
equal  the  convergence  of  momentum  transport.  Thus,  for  turbulent 
transfer  within  the  vegetative  canopy,  the  following  equation  results: 


where  the  proportionality  constant  S  can  be  further  subdivided  as 
follows: 


(14) 
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where  A  is  the  effective  aerodynamic  surface  area  of  the  vegetation 
per  unit  volume,  and  C*j)  is  the  drag  coefficient  of  the  leaf-stalk 
configuration.  Ordway  et  al  assumed  that  C'j^  is  independent  of 
Reynolds  number t  Further,  they  assumed  that  the  eddy  viscosity  was 
constant  within  the  vegetative  canopy.  In  a  later  report.  Tan  and 
Ling  allowed  for  a  variation  of  the  eddy  viscosity  and  suggested 
that  it  might  increase  linearly  with  height  as  it  does  in  the 
surface  boundary  layer.  It  is  at  this  point  where  our  views  are 
different  from  those  of  Tan  and  Ling.  Through  equation  (1)  and  (2) 
we  saw  that  K  ■  u^^  .  In  the  surface  boundary  layer,  u*  is  nearly 
constant  while  K  and  /  increase  linearly  with  height.  However, 
down  in  the  canopy,  changes  in  K  should  be  largely  due  to  changes 
in  u«  since  due  to  the  restrictive  action  of  the  leaves  and  stems 
the  mixing  length  within  the  canopy  could  be  nearly  constant.  Con¬ 
sequently,  we  set  out  to  determine  the  size  and  variation  of  the 
mixing  length  within  a  vegetative  canopy. 


In  consonance  with  the  ideal  canopy  previously  defined,  let 
us  assume  the  following  equation  to  be  valid  for  the  mean  velocity 
within  the  vegetative  canopy: 


In  U)  h  W  -  a 

dt  “‘U  az 


(15) 


where  a  is  a  constant  and  h  is  the  height  of  the  canopy.  If  x  is 
defined  as  a  dimensionless  height  equal  to  z/h  then  from  equation 
(15)  we  further  observe  that  a*  (In  U)/Sx*=0  and  U  =  Uj^exp  [a(x-l)] 
where  is  the  mean  wind  velocity  at  the  top  of  the  canopy. 
Remembering  that  i  „  ^»(ln U)  r^(ln  U)p.  we  find  that  equation 

U  ^  I  bx  I 

(13)  reduces  to: 

aK  =  hi§y  »  hiSUh  exp  [  a(x-l)  ] 

9x  a  a  (l6) 


Equation  (l6)  can  be  solved  analytically  with  the  result  that 

KU  =  KhU^j  -  hfS  dig  -  U8) 

2a* 

(17) 

Since  U*  =  K]|IJ  ^  KUa  ,  we  find  the  following  for  the  friction  veloc 
9x  h 

ity: 

“  (u*)h  -  M  (U2  -  u*) 

2.  h  (18) 


n 


Farthernore,  U*  aooordlng  to  equation  (1)  ao  /q  » 

“TT" 


(19) 


where  =  2[u*h/UjQ*  which  is  the  groee  drag  coefficient  of  the 
vegetative  canopy  ae  seen  by  the  surface  boundary  layer. 

The  definition  of  a  drag  coefficient  for  the  surface  boundary 
layer  has  proved  fruitless  in  the  past,  since  its  magnitude  depended 
very  much  upon  the  height  at  which  the  velocity  was  measured.  However, 
it  is  seen  that  the  drag  coefficient  defined  for  the  top  of  the  ideal 
canopy  has  special  significance  since  (u«/U)  is  constant  within  the 
ideal  canopy. 

In  discussing  equation  (19),  let  us  assume  that  S  is  constant 
with  height  within  the  canopy.  From  equation  (12)  we  see  that 
essentially  we  are  assuming  that  the  leaf-stalk  area  and  drag 
coefficient  functions  are  constant  with  height.  This  is  in  further 
agreement  with  the  ideal  canopy.  On  the  basis  of  similarity,  we 
shall  assume  that  the  mixing  length  is  also  constant  with  height 
within  the  canopy.  With  this  condition,  the  following  relationships 
result  from  equations  (14)  and  (19): 


Cd  2Cd  2  Cd 

(20a) 

h^s  ,  1^4  .  iiuspsc  =  ,3 

(20b) 

01 

II 

(20c) 

sic=  =  2  Cd  Va 

2  2 

(20d) 

where  LSC  ■  hA.  Our  Leaf -Stalk  Configuration  term,  LSC,  is  closely 
related  to  the  Leaf -Area-Index  (LAI)  often  quoted  ^  agronomists.  The 
LAI  is  the  total  leaf  area  of  a  ei*op  per  unit  of  ground  area.  For 
a  mature  rigid  agricultural  croo,  the  LAI  is  usually  very  near  to 

3.5. 
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Equations  (20a-20d)  will  be  considered  the  basic  relationships 
which  oharaoterise  the  "ideal  canopy".  These  relationships  could 
also  have  boon  derived  directly  from  equation  (12)  if  C’^/L  were 
tAddO  ^0  is.  To  reiterate t  the  ideal  canopy  is  character¬ 

ised  by:  (1)  there  is  uniform  vortical  distribution  of  both  the 
area  density  and  the  drag  coefficient  of  the  leaf-stalk  configura¬ 
tion;  (2)  the  drag  coefficient  of  the  leaf-stalk  configuration 
is  Independent  of  local  Reynolds  number;  (3)  the  mean  velocity 
distribution  within  the  canopy  is  exponential  and  there  is  (4)  a 
constant  mixing  length  and  (5)  a  constant  turbulence  intensity  exist¬ 
ing  within  the  canopy. 

These  characteristics  of  the  ideal  canopy  should  be  compared 
with  the  "ideal"  surface  boundary  layer  as  derived  from  the  mixing- 
length  hypothesis,  that  is,  the  logarithmic  wind  profile  e^qjressed 
by  equation  (4),  In  this  case,  the  velocity  profile  is  the  opposite 
of  that  of  the  ideal  canopy;  that  is,  instead  of  being  an  exponential 
function,  the  velocity  is  a  logarithmic  function  of  height.  The 
turbulence  intensity  (u*/U),  instead  of  being  constant,  decreases 
rapidly  with  height;  and  the  mixing  length  incre-res  as  a  linear 
function  of  height. 

In  deriving  the  equations  of  the  ideal  canopy  (20a-d),  it  was 
necessary  to  combine  the  characteristics  of  these  two  idealized 
models.  It  was  assumed  that  the  velocity  and  its  first  derivative, 
and  the  mixing  length  are  continuous  functions  across  the  plane 
defined  by  the  top  of  the  canopy.  On  the  basis  of  this  assumption, 
it  was  possible  to  define  the.  gross  canopy  drag  coefficient  Cp  which 
is  used  in  equations  (20a),  (20e),  and  (20d)  to  relate  the  canopy 
layer  to  the  boundary  layer.  One  of  the  consequences  of  this 
assuo^tion  is  that  the  second  derivative  of  the  velocity  and  the 
first  derivative  of  the  mixing  length  are  discontinuous  functions 
through  the  plane  defined  by  the  canopy  top. 

The  velocity  profiles  which  result  from  the  ideal  canopy  are 
as  follows: 


Above  the  canopy 
InA  ^  1 


Within  the  canopy 
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These  equations  are  analogous  to  the  logarithmic  wind  profile  based 
on  the  roughness  and  zero-plane  displacement  In  that  three  parameters 
are  required  to  specify  the  wind  profile  equation  (6)  ,  With  regard 
to  the  above  equations »  these  parameters  are  with  a  combination 
of  any  two  of  the  following:  /g,  S  or  Cjj* 

In  section  2,  It  was  noted  that  the  roughness  length  and  zero- 
plane  displacement  were  not  conservative  properties  of  tall  vegeta¬ 
tion.  The  profile  equations  of  the  Ideal  canopy  can  be  an  Itqprovement 
only  If  the  parameters  S,  Cd,  and  £  are  either  conservative 
properties  or  In  some  manner  unique^  Interdependent  for  a  given 
height,  density,  and  type  of  vegetation.  Equations  (20a-d)  are 
Inadequate  to  substantiate  this  condition  and  thus  recourse  to  exper¬ 
imental  data  Is  required. 


REAL  CAtJOPIBS 

In  the  previous  section,  a  hypothetical  model  for  turbulent 
transfer  within  a  vegetative  canopy  was  derived.  The  cotqparlson  of 
this  model  to  real  canc^y  flows  Is  severely  limited  by  e  paucity  of 
data.  However,  some  data  are  published  and  this  will  allow  comparison 
of  the  Ideal  canopy  with  real  canopies. 

One  of  the  basic  assumptions  of  the  Ideal  canopy  was  that  the 
average  leaf-stalk  configuration  is  essentially  invariant  with  height 
within  the  canopy.  Figure  1,  published  by  Allen,  Yocum  and  Lemon 
(Reference  9),  shows  the  accumulative  I^eaf  Area  Index  (LAI)  in  a 
mature  corn  field  indicating  considerable  uniformity  of  the  height 
distribution  of  leaf  area.  Then  as  previously  stated,  the  parameter 
A  Introduced  in  equation  (14)  is  closely  related  to  the  volumetric 
leaf  density  which  Is  the  slope  of  the  line  shown  In  Figure  1  •  As 
can  be  seen,  the  slope  Is  reasonably  constant  over  a  major  portion  of 
the  canopy  depth.  Figure  1  represents  only  the  leaf  area;  If  the 
stalks  and  adventitious  roots  at  the  base  and  the  tassles  at  the  top 
were  Included,  the  plot  of  figure  1  would  indicate  less  curvature  and 
more  closely  approximate  a  straight  line.  In  any  event.  It  Is  not 
unreasonable,  for  a  first  approximation,  to  assume  A  to  be  constant 
throughout  the  canopy  of  at  least  a  mature  com  field. 

Figure  2  depicts  wind  profiles  above  and  within  a  com  field 
canopy  as  published  by  Tan  and  Ling  (Reference  8).  Only  two  of  the 
profiles  have  data  in  the  lower  portion  of  the  canopy.  It  is  to  be 
noted  that  the  velocity  profile  is  convex  upward  within  the  canopy 
except  in  the  lowest  30cm.  The  concave  upward  profile  in  the  lower 
portion  reflects  the  dominance  of  the  ground  surface  over  the  vegeta¬ 
tion  in  controlling  the  turbulence.  Figure  3  is  a  plot  of  log  of 
wind  speed  U,  versus  relative  height  within  the  canopy  for  the  two 
complete  profiles  of  figure  2.  It  is  encouraging  to  note  that  the 
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HEIGHT 


A  MATURE  CORN  CANOPY. 


X=2/h 


data  fit  a  straight  line  reasonably  well.  Thus,  except  for  the 
region  very  near  the  ground  surface,  the  velocity  profile  within 
com,  as  i^orted  by  Tan  and  Ling,  is  exponential.  The  exponential 
profile  leads  to  a  constant  turbulence  intensity  within  the  crop. 
Observations  taken  by  Nakagawa  (Reference  10)  in  a  rice  paddy  canooy 
(see  Table  I)  also  show  the  turbulence  intensity  (u'*)2/u  nearly 
constant. 


TABLE  I 


Vertical  Distribution  of  Horizontal  Air  Flow  in  A  Rice  Paddy  - 
height  of  crop  =  90cm.  )after  Nakagawa,  Reference  10) 


Z  cm 

U  cm/ sec 

(u!ii! 

u 

40 

12.3 

0.321 

55 

13.6 

0.346 

70 

22.3 

0.331 

85 

25.7 

0.320 

100 

52.1 

0.331 

According  to  enuation  (19)  the  condition  of  constant  S  within 
the  canopy  and  an  -  xponential  velocity  profile  were  necessary  but 
not  sufficient  conditions  to  support  the  idea  that  the  mixing  length 
is  constant  within  the  ideal  canopy.  The  verification  of  constant 
mixing  length  within  the  canopy  from  field  data  is  somewhat  arduous, 
but  can  be  accomplished  by  the  following  orocedure. 


^•'e  are  concerned  with  two  distinct  regions  of  turbulent  flow  - 
that  within  the  canopy  and  that  immediately  above,  ’fe  express  the 
equations  of  momentum  transfer  in  these  two  regions  as  follows: 


K  ^  ' 

Bz 

>z 

K  ~ 

?)Z 

Bz 

=  0 

=  Su3 


Above  the  canopy 
Within  the  canopy 


These  equations  can  be  combined  into  a  single  equation  if  we  merely 
set  S  =  0  for  z>h,  where  h  is  the  height  of  the  canopy.  VJith  this 
in  mind,  we  rewrite  the  above  equations  as: 


Iz  (21) 

For  computational  purposes  it  is  necessary  to  express  equation  (21) 
in  non-dimensional  terms.  To  do  this  we  choose  the  mixing  length  and 
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velocity  at  height  2h  as  reference  values.  With  these,  we  define  the 
following  nondlmensional  variables:  x  »  s/2h,  y  =  ^/^ah  *  ^  =£Vi\h, 
Substituting  these  terms  in  equation  (21)  we  obtain: 


We  let  B 
becomes: 


=  8h3s/^%h 


8h»S..a 

r.h 


(22) 


and  further  differentiate  so  that  equation  (22) 


Equation  (23)  retains  the  unknown  parameter  B,  To  eliminate  this 
term  we  define  a  new  variable  E  =  L/B.  Using  this  variable  equation 
(23)  can  be  expressed  as  follows: 


Equation  (24)  is  the  requisite  equation  to  solve  for  the  mixing  length. 
It  must  be  remembered  that  the  variable  =  0  for  z>h. 


Equation  (24)  was  solved  for  the  two  complete  profiles  of  Figure 
2.  To  do  this  it  was  necessary  to  determine  the  coefficients  C,  and 
C, .  The  number  of  actual  data  points  contained  in  the  profiles  of 
Figure  2  is  too  few  to  adequately  evaluate  the  velocity  derivatives. 
Having  no  other  guidelines,  we  accepted  the  profile  curves  of  Tan 
and  Ling  as  the  best  fit  to  the  data.  The  portion  of  Figure  1  from 
the  ground  surface  to  2h  (556cm)  was  subdivided  by  20  equally  spaced 
grid  points  and  the  velocities  at  each  grid  point  were  extracted  from 
Figure  2,  From  those  data,  the  velocity  derivatives  were  estimated 
from  Lagrangian  polymominals.  Once  the  coefficients  C,  and  C,  were 
determined,  equation  (24)  was  integrated  from  the  ground  surface  to 
the  height  2h  using  a  fourth-order  Runge-Kutta  integration  routine* 


Figure  4  is  a  graph  of  the  nondimensional  mixing  length  (/g)  for 
the  two  solutions  of  equation  (24)  reported  by  Cionco  (Reference  11). 
In  a  general  way,  figure  4  is  compatible  with  the  ideal  canopy. 

V/ithin  the  canopy,  the  mixing  length  is  essentially  constant  except 
near  the  ground  surface;  above  the  canopy  the  mixing  length  increases 
linearly  with  height.  These  conclusions  are  not  as  well  supported  by 
the  solution  for  corn  #2  as  with  corn  #1.  However,  considering  the 
experimental  errors  which  may  be  involved,  there  is  no  reason  to 
attribute  significance  to  the  deviations  from  the  ideal  in  the  corn 
#2  solution. 


19 


RELATIVE  HEIGHT 


FIGURE  4.  THE  COMPUTED  NONDIMENSIONAL  MIXING  LENGTH  TERM  AS  A  FUNCTION 
OF  RELATIVE  HEIGHT  WITHIN  AND  ABOVE  A  MATURE  CORN  CANOPY.  (REF  11) 


Figures  3*  And  4,  In  a  qualitative  manner,  agree  with  the 
general  characteristics  of  the  ideal  canopy.  To  quantitatively 
evaluate  the  equations  (20a-d),  it  is  essential  to  compute  the 
values  of  S,  Cd  and  £  ,  This  can  be  done  by  assuming  that  the  pro¬ 
portionality  factor  between  the  mixing  length  and  the  height  is  von 
Karman's  constant.  The  nondimenslonal  mixing  lengths  were  equated 
to  linear  functions  of  height  with  least-squared  error  incorporating 
all  values  of  /e  above  the  canopy.  These  are  the  straight  lines 
shown  in  figure  4.  In  addition,  tne  average  value  of  /g  within 
the  canopy  (excluding  the  region  near  the  ground)  was  confuted.  With 
these  relationships  the  values  of  S,  the  zero-plane  displacement  D, 
and  /c  can  be  determined.  The  velocity  ration  U/Uh  is  proportional 
to  ln(Z  -  D)  by  the  factor  ;  thus,  Cp  can  be  computed 

from  a  least  squared  error  linear  function  relating  U/Uj^  to  ln(Z  -  D). 
Table  II  contains  the  computed  values  for  the  two  corn  profiles  of 
the  requisite  parameters  of  the  ideal  canopy. 

Ttm  ii  . 

_ h _ P _ a _ Cp _ 5  ic  k(h-D) 

Com  #1  278cm  257.2cm  2.000  2.462X10-*  1.875X10-“  15.83cm  8.32 
Com  #2  278cm  227.8  1.795  4.370X10-*  3.739X10-'*  27.26cm  24.08 

It  is  immediately  apparent  that  Cj),  S,  and  ig  are  not  conserv¬ 
ative  properties  of  a  mature  com  canopy.  Like  the  roughness  length, 
these  parameters  increase  with  increasing  wind  speed;  however,  it  is 
easily  shown  that  individually  their  relative  increase  is  less  than 
that  of  the  roughness.  On  the  other  hand,  the  most  conservative 
property  of  the  canopy  appears  to  be  a,  the  derivative  of  In  U  with 
respect  to  relative  height  z/h  within  the  canopy. 

TABLE  ly. 

a  hS/Cp  (h®S/24)  ^  s4  2(CD/2)»/a 

Com  fl  2.000  2.115  2.020  2. 97X10-®  2. 73X10-® 

Corn  #2  1.795  2.375  1.753  1.019X10-*  .646X10-* 

In  Table  III,  the  numerical  values  for  terms  involved  in  equations 
(20a, b,d)  are  presented.  Although  the  values  are  close  to  those 
defined  in  the  ideal  canopy,  it  is  evident  that  the  values  of  Cd  are 
consistently  too  small.  This  failure  of  the  ideal  canopy  is  due  to 
the  manner  in  which  the  ideal  canopy  was  .joined  with  the  logarithmic 
wind  orofile  of  the  surface  boundary  layer  at  the  top  of  the 
canopy. 
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In  the  above  computations  it  was  assumed  that  the  canopy  top  was 
278cm— the  figure  given  by  Tan  and  Ling.  From  a  practical  standpoint, 
it  must  be  somewhat  difficult  to  define  the  ’’top"  of  a  com  field. 

If  a  lower  reference  plane  had  been  chosen,  Cq  would  have  been 
larger,  and  consequently  the  figures  of  Table  II  would  likely  be 
more  compatible*  Nevertheless,  it  is  apparent  that  one  of  the 
difficult  problems  associated  with  the  canopy  approach  is  the  manner 
in  which  the  ideal  canopy  is  joined  to  the  logarithmic  wind  profile. 

In  developing  the  ideal  canopy,  we  did  not  consider  any  variation 
of  S  with  height,  and  it  is  at  the  top  where  this  is  likely  to  be 
most  significant.  This  is  well  illustrated  in  Figure  5  which  presents 
solutions  of  com  #1  for  various  assumed  shapes  of  the  S  function 
at  the  canopy  top.  Table  IV  lists  the  magnitudes  of  the  canopy 
parameters  which  result  from  these  profiles.  These  results  clearly 
reveal  that  the  mixing-length  profile  is  sensitive  to  changes  in 
shape  of  the  S  function  and  that  1^4  k(h  +  D). 


S 

table  IV 
/c 

’iS/Cd 

a 

Com  #1A 

2.230X10-* 

17.26 

2. 462X10- a 

2.52 

2.00 

Com  #1B 

1. 875x10-* 

15.83 

2. 462X10- a 

2.115 

2.00 

Com  #1C 

I.55OXIO-* 

14.39 

2.462X10-* 

1.75 

2.00 

The  results  for  com  have  clearly  shown  that  S  and  are  not 
consei^ative  properties  since  both  increase  with  increasing  wind 
speed.  The  increase  of  S  with  wind  speed  might  be  attributed  to 
increased  effectiveness  of  the  drag  of  the  leaves  due  to  flapping. 
However,  this  appears  to  invalidate  the  assumption  that  the  local  drag 
coefficient  of  the  vegetation  is  independent  of  Reynolds  number.  This 
does  not  invalidate  the  ideal  canopy  but  it  does  complicate  the 
equations.  Annex  A  contains  the  derivations  of  equations  analogous 
to  equations  (20a)  thru  (20d)  which  include  a  local  Reynolds  number 
effect.  They  show  that  the  mixing  length  and  turbulence  intensity 
cannot  be  constant  within  the  canopy.  Since  Tan  and  Ling's  data 
taken  in  mature  com  and  Nakagawa's  data  taken  in  a  rice  paddy  both 
show  a  constant  turbulence  intensity  within  the  canopy,  it  was 
assumed  the  Reynolds  number  effects  could  be  neglected  for  crops  of 
this  configuration,  structure,  and  elasticity.  Undoubtedly,  these 
features  are  related  to  the  maturity  of  the  two  crops  and  the 
assumption  may  prove  invalid  for  less  mature  fields. 

It  is  interesting  to  investigate  the  density  variation  of  rough¬ 
ness  elements  on  canopy  characteristics.  The  only  quantitative  data 
available  for  t'rds  purpose  is  the  bushel  basket  data  of  Kutzbach 
(Reference  3)»  Figure  6  shows  the  relationship  between  the  density 
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RELATIVE  HEIGHT  X»z/h 


MIXING  LENGTH  TERM  ,  Wl 


5.  SOLUTIONS  OF  THE  MIXING  LENGTHS  IN  MATURE  CORN  AS  A  FUNCTION 
OF  THE  SHAPE  FACTOR,  S,  FOR  CORN  . 


function  A  with  the  drag  eoeffioient  at  basket  height  less  the 
drag  coefficient  at  basket  height  associated  with  only  the  ioe 
surface.  This  correction  for  the  effect  of  the  ice  is  important 
only  at  small  values  of  A.  The  least  squared  error  analysis  results 
in  the  following  relationship: 


-  ^Dice 


=  3.904  X  10-* A 


1.243 


(25) 


Utilieing  equation  (14),  we  can  write  equation  (20d)  in  the  form 
ACp/gs  2  (Ci)/2)®/^  .  At  first  thought  it  might  seem  reasonable  to 
assume  that  the  local  drag  coefficient  (Cq)  of  a  rigid  bushel  basket 
would  be  a  constant,  but  combining  this  equation  with  equation  (25) 
leads  to  the  result  that  ^  ^g,7  if  C^were  a  constant, 

this  would  lead  to  the  ridiculous  result  that  would  vary  more  than 
an  order  of  magnitude  over  the  range  of  A  given  in  figure  6.  We  can 
conclude  that  C*p  is  not  independent  of  the  density  function  A  that 
is  the  drag  of  the  baskets  increases  with  increased  packing.  This  is 
qualitatively  in  agreement  with  the  results  indicating  that  compact 
75em  high  alfalfa  is  a  rougher  canopy  than  loose  278cm  high  com.  As 
a  consequence,  only  the  parameter  S  has  significance  with  regard  to 
the  ideal  canopy. 


DISCD8SI0K 


In  brief,  the  problem  of  the  canopy  flow  study  was  to  develop  a 
canopy  flow  model  from  basic  parameters.  The  ideid  concept  told  us 
that  the  model  is  characterised  by  the  following  properties  within 
the  canopy  for  the  steady  condition  and  no  advection: 

(1)  the  mixing  length  is  constant, 

(2)  the  mean  velocity  distribution  is  e:qponential  and, 

(3)  the  turbulence  intensity  is  constant. 

Further  development  yielded  a  set  of  four  equations  that  are  con¬ 
sidered  the  basic  relationships  that  characterise  the  ideal  canopy. 
In  conibining  these  characteristics  with  those  of  the  "Ideal  Surface 
Boundary  Layer"  of  the  mixing-length  hypothesis,  the  development 
yielded  two  wind  profile  equations  that  relate  the  canopy  to  the 


■*’‘*’The  A  is  not  Kutsbach's  A*  but  rather  is  the  reciprocal  of  the 
specific  area  defined  by  Kutsbach. 


boundary  layer  by  means  of  the  total  drag  coefficient  of  the  canopy. 
The  equations  dictate  that  the  canopy  wind  profile  is  an  exponential 
function  of  height. 

It  was  entirely  fortuitous  that  the  first  attempt  to  verify  the 
ideal  canopy  met  with  at  least  qualitative  success  utilizing  the 
data  for  com  published  by  Tan  and  Ling.  Obtaining  solutions  of 
nearly  constant  mixing  length  profiles  throughout  most  of  the  crop 
at  the  onset  of  the  study  provided  much  encouragement.  However, 
there  are  other  published  data  which  are  either  partially  or  totally 
incompatible  with  the  ideal  canopy  as  to  the  distribution  of 
velocity  within  the  canopy.  Undoubtedly  there  are  many  reasons  for 
this,  a  few  of  the  more  obvious  being; 

(1)  failure  to  meet  the  assumptions  of  uniform  distribution 
throughout  the  canopy, 

(2)  neglect  of  the  Reynolds  number  effects, 

(3)  omission  of  terns  in  the  transport  equations, 

(4)  inadequate  sanpling  in  the  observations  and, 

(5)  topographic  and  natural  surfaces  irregularities  at  or 
surrounding  observation  site. 

Figures  7a-c  present  data  from  the  literature  for  different 
vegetative  covers  and  the  double  bushel  basket  experiment  of  Kutzbach. 
The  various  vegetative  covers  represented  are  rice  paddy,  brush, 
timber  forest,  sugar  beets,  wheat,  corn,  and  citrus  orchard.  In 
general,  it  is  evident  that  "a"  of  equation  (15)  is  not  constant  in 
the  canopy  as  proposed  by  the  ideal  canopy.  However,  in  some  vegeta¬ 
tion  the  variation  is  small.  A  class  of  canopy  profiles  which  is 
entirely  incompatible  with  the  concept  that  the  momentum  transport 
is  proportional  to  the  velocity  shear  is  represented  by  those  of  Fons 
(Reference  12)  in  brush,  and  those  of  Lemon  and  Stoller  (Reference  5) 
in  immature  com  and  wheat.  These  profiles  are  characterized  by  an 
extended  region  of  the  canopy  profile  where  shear  is  negligible. 

This  class  of  canopy  profiles  implies  an  ad'litional  term  in  the 
turbulent  transfer  equation. 

The  solution  for  the  mixing  length  in  the  high  wind  speed  case 
for  the  com  data  of  Tan  and  Ling  (Reference  8)  also  might  be 
construed  to  indicate  an  additional  term.  In  the  solution,  the  ■.  r  . 
transport  of  momentum  is  accounted  for  completely  by  ^u/^z.  But  r 
this  group  of  profiles  the  velocity  shear  is  nearly  zero  throughout 
most  of  the  crop,  thus  the  mixing  length  must  vary  within  the  canopy 
to  balance  the  equation.  The  mixing  length  for  the  com  #2  solution 
indicates  a  slight  maximum  in  the  profile  within  the  crop.  This 
tendency  has  also  been  observed  in  mixing-length  analyses  made  by 
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FIGUBt  7A.  OBSERVED  WIND  PROFIIES  IN  VARIOUS  VEGETATIVE  CANOPIES. 


FIGUKE  7g.  OBSERVED  WIhD  PROFILES  IN  VARIOUS  VEGETATIVE  CANOPIES. 
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FIGURE  7C.  OBSERVED  WIND  PROFILES  Iin  VARIOUS  VEGETATIVE 


CANOPIES. 


Uchi.iiraa  (Reference  14),  This  type  profile  can  be  qualitatively 
reconciled  with  the  general  theory  of  turbulent  transport,  but  the 
requisite  data  are  nonexistent* 

Intuitively  one  would  expect  a  considerable  amount  of  turbu¬ 
lence  to  be  generated  in  the  upper  portion  of  the  canopy  as  the 
wind  encounters  the  crop*  This  would  mean  that  q®  (U***  +  v*®  +  W*^  )w’ 
would  be  fairly  large  in  the  upper  portion  of  the  canopy.  Since  this 
term  must  be  zero  at  the  surface,  the  gradient  within  the  mid  portion 
must  be  significant.  How  wire  anemometer  measurements  made  by 
Stoller  (Reference  5)  seem  to  indicate  this  is  true.  If  one  makes  an 
analogy  between  the  second  order  momentum  fluctuations  with  the 
second  order  temperature  fluctuations  as  discussed  by  Deardorff 
(Reference  1^),  it  would  appear  that  the  flux  due  to  the  diffusive 
t-erm  ^  ^  w  +  ^  pw  )  is  in^ortant  especially  at  moderate  wind 

2  p 

speeds. 

For  those  crops  that  are  especially  uniform  in  configuration,  the 
mixing  length  should  be  quite  restricted  in  its  growth  and  be  fairly 
constant  throughout  most  of  the  crop.  The  question  may  arise  as  to 
why  a  unique  mixing  length  does  not  exist  within  the  canopy,  especially 
when  it  is  emphasized  *that  the  ideal  canopy  has  a  uniform  distribution 
of  the  area  density  and  the  drag  coefficient  of  the  leaf-stalk  config¬ 
uration.  The  plant  configuration  should  also  restrict  the  size  of 
the  eddies.  If  the  vertical  area  distribution  and  the  physical 
spacing  of  the  plants  are  the  only  controlling  factors  of  the  size 
of  the  eddies,  then  the  mixing  length  within  the  canopy  should  remain 
fixed  regardless  of  the  wind  speed.  However,  the  solutions  for  the 
mixing  length  in  the  real  canopy  show  nearly  a  two  fold  increase  for 
an  increase  in  the  wind  speed  of  233cm/sec  to  317cm/sec  at  the 
canopy  top.  One  should  expect  an  increasing  mixing  length  since  the 
resulting  basic  relationships  of  the  ideal  model  indicate  that  the 
mixing  length  in  the  crop  is  a  function  of  the  coefficient  of  the 
total  drag  force  being  exerted  upon  the  canopy  volume.  The  drag  in 
turn  is  a  function  of  the  air  flow  distribution  above  the  canopy, 
which  governs  the  degree  of  penetration  of  air  flow  down  into  the 
canopy.  The  analysis  of  the  velocity  profiles  above  the  cornfield 
using  equation  6  shows  that  the  zero-plane  displacement  decreased 
with  increasing  wind  speed.  This  is  also  true  of  a  wheat  field  -  up 
to  a  critical  wind  speed.  Apparently  the  opening  of  the  canopy 
occurs  as  the  increasing  wind  speed  orients  the  leaves  along  the 
flo^M.  The  observations  indicate  that  the  mixing  length  apparently 
increases  uniformly  throughout  the  crop.  However,  at  some  greater 
speed  the  drag  exerted  on  the  canopy  reaches  a  critical  value  and  re¬ 
mains  somewhat  constant  for  greater  velocities.  In  tuini  the  size  of 
the  eddies  must  reach  a  maximum  value  similar  to  that  imposed  by  the 
plant  spacing  and  vertical  area  distribution.  Therefore,  rather  than 
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expect  a  unique  mixing  length  for  the  various  real  canopies «  one 
should  be  prepared  to  accept  a  mixing  length  In  the  crop  that  Is  a 
function  of  the  air  flow  above  the  eanop7  but  with  a  limiting  crit¬ 
ical  value. 

This  study  is  an  exploratory  atten|)t  to  understand  the  turbulent 
transfer  within  a  vegetative  canopy.  The  intent  in  developing  the 
ideal  canopy  was  heuristic.  It  is  hoped  that  this  report  will 
stimulate  thought  and  investigation  which  will  ultimately  lead  to  an 
expression  of  the  aerodynamic  roughness  effects  in  terms  of  vegeta¬ 
tive  characteristics. 


CONCLOSIOWS 

Turbulent  transfer  within  real  vegetative  canopies  is  in  general 
not  ideal.  The  utility  of  the  canopy  approach  has  not  been  proven, 
but  it  does  seem  to  show  promise.  The  relationships  developed  can 
be  applied  to  vegetative  covers  with  characteristics  that  approach 
the  assumptions  of  the  ideal  canopy  with  reasonable  success. 

The  character  of  the  mixing  length  and  the  parameter  S  with 
regard  to  the  density  and  structure  of  the  vegetation  and  its 
relationships  to  the  Reynolds  number  should  be  investigated.  It  may 
be  possible  to  do  this  in  a  wind  tunnel  or  with  controlled  e:q)eriments 
such  as  the  bushel  basket  or  in  real  vegetative  canopies.  The  manner 
in  which  the  canopy  flow  joins  the  surface  boundary  layer  flow  at  the 
top  of  the  canopy  is  another  major  problem  requiring  investigation. 
But,  foremost,  the  theoretical  approach  must  be  experimentally 
verified. 

Future  development  of  the  ideal  concept  requires  the  collection 
of  a  set  of  data  on  the  mean  velocities  and  turbulent  structure 
within  and  above  the  canopy. 
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ANNEX  A 


P.SYNOLDS  NUMBER  EFFECTS 


Previous  investigations  and  our  work  up  to  this  point  have 
assumed  that  the  flow  within  the  canopy  is  independent  of  local 
Reynolds  number.  In  certain  types  of  vegetation  or  within  certain 
limits  of  wind  speed,  the  assumption  may  be  valid.  However,  if 
one  considers  extremely  limber  vegetative  covers  such  as  alfalfa 
this  assumotion  may  be  erroneous.  The  fact  that  S  and  are  not 
conservative  properties  appears  to  invalidate  the  assumption  that 
the  local  drag  coefficient  is  independent  of  Reynolds  number.  As 
was  pointed  out,  this  does  not  invalidate  the  ideal  canopy,  but  it 
does  complicate  the  equations.  To  show  this,  let  us  define  a 
local  Reynolds  number  in  the  canopy  as  Re  =  U/Av  and  express  the 
parameter  S  as  S  =  sRe’^,  where  s  is  a  constant  and  v  is  the  kine¬ 
matic  viscosity.  Thus,  equation  (19)  becomes: 


h  r hsP.e" 
a?  (2-»-n)a 


(26) 


Equation  (26)  is  analogous  to  equation  (19) »  but  it  has  the  notable 
difference  that  it  is  not  possible  for  the  mixing  length  and  the 
turbulence  intensity  (u^/U)  to  be  constant  with  height.  If  once 
again  we  assume  a  dynamic  similarity  between  and  the  function  S, 
then  the  following  conditions • result  from  equation  (26): 


hSh  ,  h2/R2gV= 

%  CdVLSC^  2 


(27a) 


hfS, 


=  a" 


(27b) 


a 


a 


^h  J^h 


,  (2+n)  (Cj3/2)®/a 


(27c) 

(27d) 


where  Reg  is  the  gross  canopy  Reynolds  number  defined  as  Reg=(hUj^/v  ), 
Equations  (27a-d)  are  very  similar  to  equations  (20a-d);  the  most 
significant  difference  is  the  inclusion  of  th®  factor  (2+n)  in  place 
of  2.  Because  of  their  greater  generality,  we  shall  now  assume  that 
the  ideal  canopy  is  defined  by  equations  (27a-d). 
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Our  previous  analysis  was  based  upon  only  two  profiles  in  a 
com  field.  Most  available  data  deal  only  with  the  flow  above  the 
canopy  and  the  results  are  usually  published  in  tenns  of  the  friction 
velocity,  zero-plane  displacement  and  roughness  length.  It  is 
possible  to  put  this  form  of  data  into  use  if  we  assume  /h  ” 
k(h-D)  and  Cj)  =  2(k/ln  (h-D)/zQ)®  and  determine  a  functional 
relationship  between  In  and  In  Cjj  and  between  In  Cjj  and  In  Uj^. 

The  (^3+  for  com  in  Figure  2  as  published  by  Tan  and  Ling 
(Reference  d)  and  data  for  alfalfa  published  by  Stoller  and  Lemon 
(Reference  5.)  were  analyzed  in  this  manner.  Least-squared  error 
linear  functions  were  evaluated  with  the  following  results; 

i^h  =  687.2  (Cd/2)^'^^^^;  Cd/2  =  1.066  x  10-®  CORN 

^h  =  48.48  (Ci3/2)*^^®^;  Cd/2  -  1040. ALFALFA 


If  these  equations  are  incorporated  into  equation  (27d) ,  the  following 
relationships  result: 

Sj,  =  3.8658  X  10"®  (Cj3/2)*^'^^;  =  2.6638  x  10’*  Uh’'^^^C0RN 

Sh  =  5.014  X  10'®  (Cd/2)‘®®^^;  =  .9251  ALFALFA 

These  results  clearly  reveal  the  striking  difference  that  exists  in 
drag  characteristics  between  a  timber,  dense  canopy  such  as  alfalfa 
and  a  semi-rigid  open  canopy  such  as  com. 

Admittedly  the  method  of  analysis  is  crude.  However,  the  results 
should  represent  qualitative  indications  of  the  Reynolds  number  effects. 
The  n  value  for  com  is  of  special  interest  since  one  usually  expects 
the  Reynolds  number  exponent  to  be  a  negative  value.  It  is  obvious 
that  more  comprenensive  measurements  of  the  wind  flow  in  and  over 
vegetative  canopies  must  be  accumulated  before  the  Reynolds  number 
effects  can  be  expressed  in  a  quantitative  manner. 
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